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Abstract 

The concept of Feshbach resonances developed for quantum mechanical scattering is applied in 
the analysis of classical light scattering off photonic crystal slabs. It is shown that this concept 
can be realized almost perfectly in these systems. As an application guided- mode resonances in 
the grating waveguide structure (GWS) are studied in detail. Using simple resonance dominance 
approximation the characteristic properties of isolated Feshbach resonances in light scattering are 
exhibited. Formation and interaction of overlapping resonances are investigated. The relevant 
parameters of the GWS are identified which control the shape of the refiectivity of interacting 
resonances as well as the enhancement of the electromagnetic field. The differences in the properties 
of TE and TM resonances is emphasized for both isolated and interacting resonances. 
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I. INTRODUCTION 



A general property of scattering of waves, e.g., of light or matter waves is the appear- 
ance of resonances. Resonances are observed in acoustic waves as well as in scattering of 
waves associated with elementary particles at wavelengths ranging from meters to 10~^^m. 
Crudely speaking there are two classes of resonances, the shape (or potential) resonances 
and Feshbach resonances. Shape resonances occur in a variety of classical and quantum 
mechanical systems in which wavelength and the size of the resonator or target are of the 
same order of magnitude. 

The concept of Feshbach resonances [H |2] on the other hand applies to quantum me- 
chanical scattering on many-body systems such as atomic nuclei |3j and has found recently 
important applications in atomic and molecular physics [HE]- A Feshbach resonance oc- 
curs if the kinetic energy of the incident particle is close to an almost stable intermediate 
"molecular" state. In the field of cold atoms, it has been instrumental that the condition 
for appearance of Feshbach resonances can be manipulated by tuning the strength of an 
external magnetic field to which the magnetic moments of the atoms are coupled. 

In this work we will show that Feshbach resonances occur in the very different context of 
scattering of light off photonic crystal slabs. The fundamental mechanism for the formation 
of resonances is the process of turning a bound state (guided mode) into a resonance "state" 
(barely radiating mode). 

We will focus our studies on a particular photonic crystal slab - the grating waveguide 
structure (GWS). Guided mode resonances in GWS have been studied both theoretically 
and experimentally since 1985, Refs. [HI |7j, and have already been incorporated into a wide 
variety of applications, cf. Refs. [814T3]. Theoretically they have been investigated using sev- 
eral different approaches including a scattering matrix approach [H] , temporal coupled mode 
theory [l5l[T6], Wigner-Weisskopf formalism [IT] as well as coupled mode theory [lIS]. In the 
majority of these theoretical approaches parameter fitting is required to make comparison 
with exact numerical calculations fl9\ 120] . Our approach is similar to the coupled mode 
approach of Rosenblatt et al. [TH] which was in turn based on the work of Kazarinov et al. 
in distributed feedback lasers [21] ■ Yet the origin of our treatment is different being based 
on the formalism of Feshbach resonances. In this way are able to derive accurate analytical 
expressions for essentially all resonance properties. Miroshnichenko [22l [23] and others [21] 
have pointed out that various photonic resonances can be treated as Fano-Feshbach reso- 
nances, yet since no analytic expressions for the coupling exist a parameter fitting procedure 
had to be be used. In addition most of the theories concentrated on the properties of the 
TE resonances. To the best of our knowledge non of the fully analytic theories of guided 
mode resonances treated the TM polarized resonances and their differences from the TE 
resonances. 

As in the case of atomic physics, the light scattering off GWS will be shown to be tun- 
able in a variety of ways generating isolated as well as overlapping Feshbach resonances. We 
will apply a simple resonance dominance approximation and will obtain essentially analytic 
expressions for all relevant observables. In this way we will provide the tools to produce 
resonances with properties desired in applications by adjusting the the GWS parameters. 
The resonance dominance approximation will also allow us to study systematically the in- 
teraction of two or more Feshbach resonances. 
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II. RESONANCE SCATTERING OF LIGHT - FORMAL DEVELOPMENT 



In this section we will develop a formalism which will make explicit the role of Feshbach 
resonances in light scattering off photonic crystals. For this purpose and as an example we 
will consider the so called grating waveguide structure [IB] (GWS) which consists of a plane 
waveguide with one dimensional grating on the top of it. The left part of Fig. [l] shows an 
example of a GWS with piecewise constant grating layer where the waveguide, the grating, 
the substrate and the superstrate layers are shown as well as the incident, reflected and 
transmitted light. The right part of the figure shows the corresponding dielectric function 
where we replaced the grating layer e// {x) with a homogeneous layer with an effective 
dielectric constant eg. The exact definition of eo depends on the polarization of the incident 
light, cf. Eq. ^ and Eq. (33) below. In our formal studies we will treat on equal footing 
both piecewise constant as well as continuous transitions of e between the dielectric layers. 

In the present work we investigate only the case of classical incidence where the incident 
light wave vector is perpendicular to the grating grooves, i.e. ky = . We start with the TE 



polarization. The modifications needed for the TM case will be considered in Section II B 
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FIG. 1: Left: Light incident on a grating waveguide structure (GWS) composed of piecewise 
z— independent layers with dielectric constants ej-...e/y-The geometry of the grating layer (II) is 
characterized by the grating period A and the duty cycle d/A. Right: Polarization dependent 
effective dielectric constant (cf. Eqs. ([6]) and (33)) as a function of z for the GWS on the left. 



A. TE waves 



In the TE polarization the electric field of the incident light is parallel to the grating 
grooves which in our geometry (cf. Figjl]) implies 

E{r) = E{x,z)ey. (1) 

The stationary Maxwell equation for a time harmonic electric field of frequency u (with the 
speed of light set to unity) 

V X V X E{r) = e{r)uj^E{r) , (2) 

then simplifies to 
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To describe light incident on GWS from z — )■ — cxo with the wave vector in the x-z plane 
k = (fc^, 0, kz) we impose the boundary conditions that for z — )■ oo we have only an outgoing 
wave. With this choice the asymptotic behavior of E[x, z) is 

lim E{x, z) = e*'=-^(e'''" + r{u, k^)e~^^^ , lim E{x, z) = t{uj, k.^^e'^'"'' e^^^ . (4) 

z— >— oo z—>ca 

The choice of the normalization is free and, as is common, we have normalized the amplitude 
of the incident wave to be 1. As a consequence, the reflected and transmitted fields are 
given in units of the incident field. The quantities r{u, k^) and t{uj, k^) are respectively the 
reflection and transmission amplitudes. One of our goals is to devise a simple method to 
calculate these amplitudes and to show that their resonant behavior is typical of Feshbach 
resonances. Note that k1 + {k~Y = e\UJ^ while A;^ + (fc^)^ = eiycu^ 

The x-dependence of the dielectric constant is limited to the grating layer {z\ < z < Z2). 
It is periodic with the grating period A 

e(a; + A, 2;) = e(x, . (5) 

and can be represented in terms of its Fourier-components 



1 



A/2 



e^{z) = - / c/xe(x,z)e-*"^«". (6) 

^ J -A/2 



Here we consider real valued e(x,z 



e.^iz) = t^iz) . (7) 

If e{x, z) is symmetric or antisymmetric around an appropriately chosen center of the ele- 
mentary interval, the Fourier-components e„(z) are real or imaginary respectively. In the 
coordinate system shown in Figjl| the periodicity of e{x,z) implies that the x dependence 
of E{x, z) contains only Fourier components with the discrete x-components of the wave 
vectors of the form 

27r 

k^ + nKg, Kg=—. (8) 

Inserting 

E{x,z)= E„(z)e-*('=-+"^«)% (9) 

n=— 00 

into Eq. ([s]) we rewrite the latter as a system of ordinary differential equations (here and in 
the following we will suppress the summation limits ±00) 

( - ^ + (A;.. + nKgf - eo{z)uj^)E„{z) = uj^ Yl ^n-Uz)Em{z) . (10) 

To see what the asymptotic conditions ^ imply for this system we note that for m ^ n 
en-m{z) = when z — )■ ±00 so that the equations decouple at large \z\. We will further 
assume that we deal with the so called subwavelength zero order grating, i.e. such that 

- (k^ + nKgf > only for n = 0. 
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Therefore all the components En{z) with n 7^ are required to decay exponentially with 
z — )■ ±00. The non vanishing asymptotic conditions Q apply only to Eq{z) 

lim Eo{z) = e'''^' + r{u, K)e-'^^' , lim ^o(^) = t{uj, k^)e'^^' . (11) 



We convert (10) into a system of integral equations and introduce to this end the Green's 



functions satisfying 

^2 + {K + nKgf - eo(2:)a;') (?n(^, z') = -6{z - z') , (12) 



dz^ 



and obtain 



E„(z) = 4+Hz)(5„o - u;' y" dz'g,,{z, z') ^ e„_^(/)E„(/) . (13) 



Here we use the fact that light is incident only in the n = channel. The term E\^\z)5no 



n - 



satisfies (10) with vanishing right hand side so that Eq \z) is a solution of 

[-^,+kl-e,{z)u')E^^\z)=Q. (14) 

and we impose the same boundary condition as in Eq. (|4]) with reflection and transmission 
amplitudes which we will denote by ro(a;, k^) and toiuj, k^) respectively. It is the z-dependent 
part of the electric field i?Q'''''(z)e*'^"'^ which propagates in and is scattered off the effective 
dielectric structure defined by ^q{z) shown in the right part of Fig.jlj We will refer to 
this scattering as "background scattering" upon which Feshbach resonances are formed by 
coupling to guided modes. In Appendix VI B the appropriate "background" Green's function 
is explicitly constructed. 

Systems of equations of coupled channels such as Eqs. (10 , 13) are underlying the descrip- 
tion of Feshbach resonances in atomic and nuclear physics. A Feshbach resonance occurs 
if, by neglecting certain couplings, a bound state of the entire system exists. Accounting 
for the couplings converts this state into a resonance. In this spirit let us consider the above 
system (10) in absence of the coupling e„(z) =0, n 7^ of the Fourier components En{z). In 
this approximation, the photonic crystal slab is described by the effective dielectric constant 
in the right part of Fig.[T] and the resulting guided modes are the progenitors of Feshbach 
resonances. 

As the wave functions of bound states in quantum mechanics, guided modes in dielectric 
slabs are localized in the z-direction 

z ^ ±00 , E{x,z)^0, (15) 

and are determined by the eigenvalue equation 

/ \ 

with the frequency dependent eigenvalues r] = rj^u). The electric field is given by 

E{x,z) =e"^''£^{z). (17) 
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Inserting this ansatz into the wave equation (cf. Eq. (|3 

^2 52 



0. 



(18) 



the dispersion relation between frequency u and the x-component /3 of the wave vector is 
obtained 

/S' = -^(w) • (19) 

In Fig. ([2]) we show as an example the dispersion curves of the first three TE and TM 
guided modes for a medium with the z-dependent dielectric constant 60(2;) of Fig. ([T]) and 
numerical values of the dielectric constant given in Eqs. (57) and (58) 
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FIG. 2: Guided mode frequencies iOi{f3) as a function of the propagation constant f3 for a medium 
with the z-dependent dielectric constant eo{z) of Fig. ([T]). The numerical values of the dielectric 
constants are given in Eqs. (57) and (58) and i is the mode order. Dashed curves are for TE 



guided modes (red color online), solid curves TM (blue color online). Also shown are the lightlines 
w = /S/y/e for e/, em and ejv- 



The subset of equations (|13) with n ^ can be rewritten as 



EJz) = - 



-Vi - [K + nKgY 



+ Gc{z, z) \^^^^ en- 



(20) 



where the first term in the curly brackets is the contribution of the guided modes S^^^ (z) to the 
the Green's function gn{z,z') (cf. Eq. (12)) and Gc{z,z') the contribution of the continuum 
(radiating) modes. The guided modes are described by normalizable functions and we require 



dz£^Xz)^*{z) 



(21) 



With this exact reformulation of the original wave equation we are in the position to formu- 
late the criterion for appearance and dominance of Feshbach resonances. Theoretically, an 
isolated Feshbach resonance is expected to occur for sufficiently small coupling {en{z) n 7^ 0) 
and if the kinematics {kr^ and u) is chosen such that (cf. Eq. (20)) 



rii + + nKgf ^ . 



(22) 
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In this case the scattering of hght will be dominated by one of the terms in the discrete 
part of the spectrum, i.e., by one of the guided modes localized in z. The coupling to the 
extended mode delocalizes these modes and turns them thereby into resonances, i.e., guided 
modes become "slightly radiating". Narrow subwavelength resonances in this system were 
observed (cf. Refs. [TT l [T8 | I251 - I27] ) in a wide range of the parameters, i.e. wavelength, angle 
of incidence and the dielectric structure parameters. This supports the assumption of weak 
coupling of the guided to the extended modes. 



The criterion ( 22 ) is essentially identical to the quantum mechanical criterion for Feshbach 
resonances which requires that the incident energy coincides approximately with that of a 
"molecular" state which in the limit of vanishing coupling becomes a true bound state. As in 
scattering of ultracold atoms, the properties of the Feshbach resonances can be manipulated 
in light scattering as well. They can be tailored towards a particular application by variation 
of the parameters of the system or externally by using electro-optical materials |28]. 

For the analysis of the resonances observed in scattering of light off photonic crystal slabs 
the reformulation of the wave equation (|3| by the system of equations (([l2|), (20)) suggests 
to neglect for tt, 7^ the continuum contributions to the Greens functions 5f„(z, z'), i.e., to 



approximate Eq. ( 20 ) by 



EM - / dz' J2 + + E ^n-M)EM) , n^O. (23) 



We will refer to this approximation as "resonance dominance". In general the condition (22) 
is satisfied for one guided mode only and the contributions from the other guided modes 
can be neglected. With a a judicious choice of the parameters, a simultaneous excitation of 
two or more Feshbach resonances can be achieved and the interaction of these "overlapping" 
resonances can be studied. A particular class of two interacting Feshbach resonances will be 
discussed later. 

Notwithstanding the formal connection between Feshbach resonances in quantum me- 
chanical scattering of particles and in classical scattering of TE polarized light the content 
of the corresponding wave equations is very different. We mention in particular the role of 
the frequency. The Maxwell equation ^ implies that the frequency plays a twofold role. On 
the one hand, both in quantum mechanics and in electrodynamics, the frequency determines 
the asymptotic properties of the incident particle or light. Simultaneously in electrodynam- 
ics, the frequency also determines the strength of the interaction of light with the dielectric 
medium. Peculiar consequences are the frequency dependence of the eigenvalues rj and the 



electric fields of the guided modes £r){z) (cf. Eq.(16)) as well as the independence of these 
quantities on the x-component of the wave vector. As we will see below, these properties 
not only distinguishes TE modes from quantum mechanical waves but also from the TM 
modes. 



B. TM waves 

In many aspects the treatment of the TM polarization is identical to the TE case so we 
will present it briefly. In the TM polarized waves the magnetic field is parallel to the grating 
grooves which in our geometry means that 

U{r) = H{x,z)ey. (24) 
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Inserting this into the Maxwell equation (we set c=l) 



V X 



-V X H = a;^H 



e r 



we obtain 



e[x,z) 



H(x, z) = uj H(x, z) . 



(25) 



(26) 



As discussed in |29] this equation is analogous to a Schrodinger equation in two dimensions 
with no potential but coordinate dependent mass. Accordingly we write it in the form 



where 



e 



(27) 



(28) 



'e{x,z) e{x,z) 

As in the TE case we need to solve the above equation under the condition that only outgoing 
wave is present at z — )■ oo 

z—^oo 



lim H{x, z) = e'^-^e"'^'' + r{u, K)e-'^^') , lim H{x, z) = t{u, h)e''-''e"'^ ' (29) 



where as before r(a;, kr^ and tiio^ k^) are respectively the reflection and transmission ampli- 
tudes for the present TM scattering. Using the analog of the expansion in Eq. ^ 



H(x.z) 



z e 



-i{kx+nKg)x 



(30) 



and inserting in (27) we obtain a set of coupled equations similar to (|10|) 

•2 Q„jHJz) 



UJ 



Z 



where we introduced 

Onrr 

and defined 



{z)dz + {kx + nKg){kx + mKg)'^n-m{z) 



ln[.Z) 



A/2 



dx- 



_^-inKgX 



(31) 



(32) 



(33) 



A J -A/2 e(a;, z) 

The effective dielectric constant ^^{z) shown in Fig.[l] is (for the TM polarization) just 

iho{z). 

We note that there exist an ambiguity in our definition of 'jn-m{z) in (32). If we view it 
as a matrix 'jnmiz:) we can either use the above definition •jnmiz) = 'Jn-miz) as given by Eq. 
(33) or we could also have defined it as the inverse of the matrix e„_m with e„ given by Eq. 
(6^ We will comment upon and test this ambiguity below in our numerical examples. 

As in the TE case the subwavelength grating condition imply that all the components 
Hn{z) with n 7^ should have exponential decay as the boundary conditions at 2; — )■ ±00. 
The zeroth component Hq{z) should asymptotically behave as 



lim Hq(z) = e 

2— 00 



+ r{u, k^)e 



lim Hq(z) = t{uj, kx)e 



ikj z 



(34) 
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Let us now convert the set (31) into integral equations similar to what we did in Eq. (13) 



H,,{z) = H^^\z)5n, + j dz'g^iz, z') ^ Qn^HM) , (35) 

where in analogy to Eq. (14| the component Hq^\z) and the Green's functions gn{z, z') solve 
respectively the differential equations 

[u'-Qoo]Ht\z) = 0, (36) 

and 

[Uj'-Qnn]gn{z,z') = 6{z-Z'). (37) 

Here Hq'^\z) describes the TM scattering off the effective structure defined by e{z) = 7(^^(2;). 



As in the TE case we choose it to satisfy the boundary conditions similar to (34) denoting the 
corresponding reflection and transmission amphtudes by ro((X', k^) and to{u, k^) respectively. 
Accordingly we chose the boundary conditions for the Green's function gQ{z,z') such that 



the full solution Hq{z) satisfies Eq. (34). 



It is useful to consider the eigenfunctions of the operators 9„n with different n's 

{z) = r]n {z) . (38) 

Note that for each n there is a complete set of eigenfunctions 'Hn,v{z) with the corresponding 
eigenvalues rjn^y distinguished by the index v. 



From the comparison with the Maxwell equation (26) one can easily see that the eigen- 
functions 'Hn,u{z) determine the z dependent part of the TM photonic modes with propaga- 
tion constant kx + nKg and frequency ^ = ?7n,i 



i^(fc.+nK,),.(x, z) = e*('=^+"^«)^-H„,,(z) (39) 

of the effective structure defined by e{z) = 7^^(z). We will call such a structure "unper- 
turbed" . We will base the treatment of the effects of the x-dependence of the grating upon 
solutions for this structure which we assume known. 



We now concentrate on the equations (35) with n ^ and use the spectral representation 



of gn{z,z') as in Eq. (20) 

J In ^ VriMi I , 

where the first term in the curly brackets is the contribution of the guided modes (denoted 
by discrete index z/j) to the the Green's function gn{z, z') while Gc{z^ z') is the contribution 
of the continuum. 

We have assumed the normalization of the guided modes 

dz'Hn,uAz)K,,Sz) = l. (41) 

Note that since 6„m with n 7^ m is vanishing outside the grating interval Ig , the integrals 



in (40) are over /„ 
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As in the TE case our main approximation will consist in neglecting the continuum 



contributions to the Greens functions 5'„(-2, z') for n 0. Thus we approximate Eq. (40) as 



Hn{z) = [ dz' ^ Q^^HUz') . (42) 

Jin ^ Vn.Ui I , 

This is the "resonance dominance" approximation in the TM case. An isolated TM Feshbach 
resonance occurs when the physical parameters are such that 

Tln,u^ ~ (43) 

and the contributions from the other guided modes can be neglected. Of special interest to 
us will also be cases of a simultaneous excitation and interaction of two or more Feshbach 
resonances. 



III. ISOLATED RESONANCES 
A. TE resonances 

In this section we will discuss in detail properties of an isolated Feshbach resonance. By 
comparison with the results of exact numerical evaluations, we will determine the degree of 
validity of resonance dominance in light scattering. To be concrete, we will carry out this 
study in the grating waveguide structure (GWS) as shown on the left in Fig.[l] The grating, 
periodic in x (cf. Eq. ([s])), is restricted in the z-direction to the region II which we denote as 
the grating interval Ig. The effective dielectric constant is given by (cf. Eq. ([6]) and the right 
part of Fig. IT]) 



eo(2:) = eiO^zi - z) + e^O^z - 2:1)^(^2 - z) + emO^z - ^2)^(2:3 - z) + ei^O^z - z^) . (44) 
while 

e„(z) = e„ e{z - Zi)9{z2 - z) for n 7^ . (45) 

In these expressions e„ denote constants given by Eq. ([6]) for the layered structure when z 
is restricted to the grating interval. 

Let us assume that these parameters of the GWS as well as the angle of illumination 



and the frequency are such that the condition (22) is satisfied for a single value oi n = v 
and a particular eigenvalue r^o of the guided mode equation (16). Let us denote by the 



corresponding eigenfunction. Beyond the resonance dominancefcf. Eq. ([23])) we truncate the 



system of equations (([T3|), (23)) further and take into account only two modes, the extended 



system of equations (13) to 



mode of the incident light (n=0) and the guided mode (r^o, n = u ). We thus reduce the 
^ to 

Eo{z) = eI;\z) - uj^ [ dz'goiz, z')E,{z') , (46) 

J la 

?7o + {K + vK, 



Eu{z) = , , £r,S^) / dz'£,,{^')E,{z') , (47) 
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and for simplicity will refer to this combined approximation to "resonance dominance". 
This system can be solved analytically. As Eq. (47) shows the resonance component of the 
electric field Ey is proportional to the resonating guided mode 

E,{z) = a,£^,{z) (48) 

with the proportionality coefficient a^, measuring the degree of the excitation of the guided 
mode. Inserting this expression for Ey{z) into Eq. (46), multiplying the resulting equation 
with Srjolz) and integrating we obtain the following expression for the field enhancement 
coefficient 



P 



(49) 



with C*^^-* denoting the coupling to the guided mode of the electric fields Eq'^\z) incident 
from either — oo or oo (cf. Eq. (153) for their precise definitions) , 



dzS,,{z)Ei,^\z) . 



The integral is carried out over the grating interval Ig. The denominator in (49) 

P = 'no + iK + ^Kgf + le^pw^s , 



(50) 



(51) 



contains the resonance condition (cf. Eq. (22)) modified by the "self-coupling" S of the guided 
mode 



dz / dz'£no{z)go{z,z')£noiz') . 



(52) 

S is generated by transitions from the guided to the extended mode, the propagation in 
the extended mode and then the back transition to the guided mode. The strength of this 
contribution to p is determined by the corresponding Fourier coefficient e^y of the dielectric 
function (cf. Eqs. (|6|, ([T]), ([45|). S is complex with the imaginary part accounting for the loss 
of intensity from the guided to the extended mode. It gives rise to a shift of the resonance 
position and to a width. Using the identity (160) derived in Appendix B, the following 
expression for the width of the resonance is obtained 

,4 



r 

2 



e,,Pa;^ ImS 



2k- 



(53) 



with the z-component of the wave vectors fcj defined inEq. (152). The strength ay of the 
guided mode excitation also determines the resonance behavior of reflection and transmission 
amplitudes which are found from asymptotics of Eq{z) in Eq. (46). Using the expression for 
(7o(-2, z') and its asymptotics derived in the Appendix B (Eq. (157)) as well Eq. (48) we find 



2k: 



i 






2k- -r]o - {k, + pKgY - 







(54) 



Similarly the corresponding expression for the transmission amplitude is obtained 
(cf. Eq. (pip 

t{u, k^) - to{uJ, k^) = ■jrz—Oy t-u C^'^ = -WTT n~, 1 — • (^^) 



2k: 
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Reflection and transmission amplitudes are related to each other. The resonance dominance 
approximation shares with the exact system of equations the property of power conservation. 
It is straightforward to derive the following identity 



2 2 

\r{uj,K)\ + -^\t{uj,h)\ = 1 



(56) 



which is valid irrespective of the number of guided modes. It thus applies to the case of 
isolated as well as overlapping resonances which we study below and is easily generalized to 
the case when more than one extended mode is present. 

The background scattering is not only present in the first term in Eq. (46) and in the 
reflection and transmission amplitudes 
and S 



To, to). It also affects the resonance terms via C*-^-* 
u, Eqs. (|50l) and (|52l). Only in this way it is possible to have the absolute squares of 



the total reflection and transmission amplitudes r{u, kx)j t{u, kx) satisfying the identity (56) 
(cf. also Eq. ( 152[ )). The interference between background and resonance contributions will 
be one source of the asymmetry of the shape of the resonance curve akin to the phenomenon 
of Fano resonances in quantum mechanical systems [301 El] • 

Fan et al. [inj have previously shown that the reflectivity curves for guided mode res- 
onance can be described analytically as Fano resonances where the spectral bandwidth of 
the resonance is an external fltting parameter from exact numerical simulations. Our inves- 
tigations support their claim and in addition produce the full analytic form of the resonant 
width (cf. Eqs. 49 55), including the weak spectral dependence of S. 

In general, the background scattering also exhibits resonances. These shape resonances 
owe their existence to the properties of the "potential" u'^eo{z) (Eq. (44)) and do not involve 
transformation of a localized mode into a (delocalized) resonance. 

We have restricted our numerical studies to GWS with piecewise constant dielectric struc- 
tures. Minor changes only are required to account for a general z-dependence. We have to 
replace 

eX^^'^ ^ j dz £no{z)ey{z) e''q^\z) , |e^|^S j j d'Z'£r,o{z)eu{z) go{z, z') e-i,iz)£r^^^{z') . 
In all our numerical studies we kept fixed the thickness of the layers and the dielectric 



constant of the superstrate. Using the notation of Eqs. (44) and (45) 



(57) 

In the first application we choose the following profile of the grating interval, the dielectric 

1 



Z2 — zi = 0.1 /im, £ = Z3 — Z2 = 0.4 /im , ei = 1. 

1 we choose the following profile of the grating intei 
constants, the grating period A, the duty cycle d/A and the grating contrast eg 

/\ [-,./ -N^/rO 0\\n/ \ n / „ . \ 

euix 



l + {eg- l)e{d^ - x^) e{x - xo)e{A + Xo-x) 

A 



em 



4, eiv = 2.25, A = 0.81 fxm, d 



eg = 4. 



(58) 



These parameters have been chosen such that, for the wavelength A = 1.5 fim and the angle 
^ = 5° of the incident light, the extended mode resonates with the guided mode with z/ = — 1 



(cf.Eqs.(46), m). 



The numerical evaluation of the expression for the refiection amplitude (54) proceeds in 
two steps. In the first step, the extended E^{z) and guided £no{^) modes together with the 
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guided mode eigenvalue rj^ are calculated for the effective medium (44) where the dielectric 
constant in the grating layer is replaced by its averaged value eo 
grating profile (58) and the parameters (58). The functions E^{; 



2.5 for the choice of the 
and Sr^Q have an analytic 



form and only the eigenvalue r]Q requires numerical solution of a transcendental equation. 
Given these building blocks it is straightforward to calculate in the second step observables 
such as the reflectivity, the electric field, etc. . 

The eigenvalue rjo as a function of u"^ and the electric fields of the guided modes are 
displayed in Fig.lsj The numerical results of the evaluation of Eq. (54) are presented on the 





0.0 , , 0.5 

z [/iinj 



FIG. 3: Results for GWS with parameters specified in Eqs.(57) 
(cf. Eq. (16)) as a function of uj'^ 



(58). 



Left: the eigenvalue r] 
Right: normalized electric fields of guided modes for A = 1.5 /im 



(solid) and 2.1 /xm (dashed line). 



left of Fig.lH 

In order to assess the accuracy of our theory we have compared it with numerical results 
obtained by applying transfer matrix techniques (cf. [3^) to the system of equations (10) 
truncated to a finite number of channels. Convergence was typically achieved with 5-10 
channels. We shall refer to these converged results as "exact". Such exact values of the 
resonance position Ares and the full width at half maximum F are 



Xr 



1.5000/im, F = 3.91168nm. 



(59) 



The deviations from these values obtained when truncating the coupled equations (|10|) to 
just two (i.e. n ■ 



-1 and n = 0) and in the resonance approximation are respectively 



SXr 



-0.6 nm, ST 



-0.0087 nm, SXr 



-0.6 nm, (5F 



-0.0129 nm. 



(60) 



The exact resonance position is shifted by 6Xres = 17.4 nm relative to the eigenvalue of the 
guided mode. The resonance dominance approximation reproduces this shift as well as the 
resonance width with an accuracy of 3.3 %. 

The coupling of the u = —1 guided mode to other guided modes neglected in this simplest 
version of the resonance dominance approximation is the most likely mechanism to account 
for the few % deviation of the resonance position from the exact result. It is not difficult to 
show that, independent of any details, the resulting shift due to such additional couplings is 
always towards longer wavelengths. In addition, the essentially identical results obtained in 
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resonance dominance and at the lowest level truncation rule out a significant contribution 
from the u = —1 continuum contributions (cf. Eq. (20)). 



Similar results and the same level of agreement have been found for other values of the 
angle of the incident light. In order to enhance the effect of the background scattering and 
to study in detail the changing shapes generated by the interference between resonance and 
background contributions (Fano resonances) we have changed to the following parameters 
of the dielectric medium (cf. Eq. (58)), 

, A = l fim. (61) 



eiv 



With this choice one finds a strong background reflection with |ro(a;, , kx) \ = 0.9. As shown 
in the right part of Fig. [4], as a result, a drastic change in the shape of the resonance 
curve is obtained due the interference between background and resonance amplitudes. To 
facilitate the comparison we have shifted the resonance dominance curve by 0.75 nm. An 
almost perfect agreement is observed indicating that the interference between background 
and resonance amplitudes is treated correctly. 





1.50 

X [|J.m] 



2.490 2.495 

X [\im] 



FIG. 4: Reflectivity (cf Eq. (54)) as a function of the wavelength. Left: angle of incidence 5° 
and parameters of the GWS specified in Eqs. (57)-(58). Solid lines: exact values, dashed hnes: 
values obtained in the resonance dominance approximation. Dotted line (in the left figure and 
coinciding with the dashed line): values obtained by truncation of Eq. (20) to n,m = 0,-1, see 



the text. Right: angle of incidence 75° with changes in the parameters (cf. Eq. (61)). Solid line: 



exact values, dashed lines: values obtained in the resonance dominance approximation and shifted 
by 0.75 nm towards larger wavelengths to facilitate the comparison. 



A peculiar feature of the TE case in sharp contrast to quantum mechanical scattering is 
the dependence of the strength of the interaction of light with the dielectric on the frequency 
~ {e{x,z) — l)a;^ (cf. E q. (|3| )). This implies that the eigenvalues and eigenfunctions of the 
guided mode equation (16) depend on the frequency (cf. Fig.js]). At the same time, as in 
quantum mechanics, the frequency determines asymptotically the behavior of the incident 
or scattered light. As a consequence of this twofold role of the frequency, the observables like 
width and position of the resonance have a rather intricate dependence on the frequency. To 
illustrate the consequences we consider the width and expand the resonance denominator p 



(51) around Ur, the zero of its real part 



Rep((X'r) = 0, p{ijj) ^ j^u — Ur — i'y ^\e^\'^ijj^lmT,{ijJr)) 



with 



7 = (Kg — ojr sin 9) sin 9 + 



2uj 



drio 



du 



fw^Re SI 



(62) 
(63) 
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In the range of frequencies of interest to us, the dominant contribution to 7 arises from the 
second term. Without a frequency dependent interaction, 7 ^ —2ur. Instead, as can be 
read off from Fig.|3] , we find 7 ~ —7ur- Thus the frequency dependence of the interaction 
leads to a significant narrowing of the resonance structure. 



B. TM resonances 



We consider here the simplest case when only one resonating term, i.e. a term with 



small denominator is present in the sums over the modes in Eq. (40). We stress that our 



assumption means that such a term occurs only for one particular value of n, i.e. that the 



resonant condition (43) is fulfilled for one particular pair (n, Ui). We denote it by {k, k). We 
approximate by truncating the system (35) to only two equations - that of = (in which 



the incoming light appears) and the resonating n = k. This is the truncation approximation. 
It leads to two coupled equations 



Ho{z) = H^o^\z)+ [ dz'go{z,z')eokHk{z' 

Jig 



QkMz') 



(64) 



(65) 



Furthermore in the equation for Hk we retain only the resonating mode fc, k. This is the 
approximation of the resonance dominance. We thus obtain 



^o(^) 



Hi^\z) 



dz'go{z,z')QokHk{z') 



Hkiz) = 2^h^ f dz'nU^') QkoHoiz') . 



(66) 
(67) 



Exactly as in the TE case this system of equations can be solved analytically, cf., Eqs. 
( 46p7 ). Equation (67) shows that the Hk{z) is just proportional to the guiding mode 

'Hk^niz) 

Hk{z) = aUkA^) 5 o- = — / dz''Hk,K{z')QkoHo{z') . (68) 



- r]k,K J I, 



Inserting Hk{z) = al-ik^n{z) in (66) we obtain 



Hq{z) = H^~^\z) + dz'go{z,z')Qok'Hk,K{^') ■ 



(69) 



We then multiply both sides of this equality by 'Hk,K{z)Qko, integrate over z and solve for 
the field enhancement coefficient a. We obtain 



a 



where we have introduced 

c(±) = / dznkMekoH!t\z) 



S = / dzdz'Hk,K{^)&ko9o{z,z')&ok'HkA^') ■ 

'la 



(70) 

(71) 
(72) 
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Here Hq ^ is the matching pair of the solution Hq^^ as defined in the TM part of Appendix 

lYlB 



Equations (68) through (72) provide a complete solution of the problem. It is expressed 
through the solutions Hq'^\z) and T-Lk,K of the unperturbed structure. From (69) one can 



easily extract the refiection and transmission amplitudes from the asymptotic behavior of 



Hq{z) and g^^z^z'). Using Appendix VI B we obtain 



lim Hq{z 



lim Hq{z) 



|2(3{+)2 



t e 



t = tn- 



2k- u;2 _ _ s ' 

2k+ uj^ - r]k,^ - S 
r(+)/ 



(73) 



Here ro and to are the "background" values which appear in Hq (z). N ote also that we 
have used the replacement 6„m — ^ In-m^nm explained in Appendix VI C where 

Qnm = X^z + {K + nKg) {k, + mKg) . (74) 

to define C^^) = -fkC^^^ with 



dz 



k,K\.z)d:,H'f-\z) + {k^ + kKg)k^Hk,n{ 



)H, 



(75) 



As in the TE case the self coupling S which appears in the denominator of (70) and (73) 



plays a crucial role. Its real part shifts the resonance frequency away from the eigenvalue 
rjkK of the resonating guided mode while the resonance width is given by the imaginary part 
of E. 

Although formally the expressions for r and t are very similar to those for the TE case 
one should note the important and crucial differences in the expressions ( |7l||72| for C(±) and 
S which contain the coupling operators Qmn- 

As already mentioned the TM formalism has an ambiguity of using jrn-niz) as given by 
Eq. (33). We could alternatively regard em-n{z) as a matrix e^n = ^m-n{z) and inverting 



it to find 7m-n = (e~ )mn- We have found that the second alternative leads to improved 
comparison with exact results and adopted it in our examples below. We do not have a good 
argument to justify this improvement. In the literature, cf., Refs [331437] . this ambiguity in 
the truncation has been discussed and resolved in favor of the "inverse" method by consid- 



ering the asymptotic convergence properties of the Fourier components of H (cf. Eq. (30)) 
which is an important issue in high precision numerical studies. 



In our examples below we have used the same parameters of the structure (57) and 
(58) as in the TE case. In the first step, the extended Hq{z) and guided 'Hk,K{z) modes 



together with the guided mode eigenvalue rjk^K were calculated for the effective medium (44). 



They were computed analytically apart from finding eigenvalue rjk^K which required solving 
numerically a transcendental equation. 

The main difference at this step as compared with the TE case is the peculiar bound- 
ary conditions for the derivatives of the solutions at the boundaries. The combinations 
^oiz)dzH^{z) and eo{z)dzT-Lk,Kiz) rather than the derivatives themselves must be continuous. 
This leads to a subtlety when the integrals (71) and (72) for S and C*-^^ are computed. Since 
the derivative terms in the operators Qko and 6ofc contain dz^k{z)dz rather than dzeo{z)dz 
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combinations they produce 5-function like singularities when acting on Hq (z), 'Hk,K and 
gQ{z,z'). We refer to Appendix VIC in which we explain how we dealt with such singu- 
larities. Adopting this approach, the quantities S and C-^-* essentially can be evaluated 
analytically with the results used to calculate reflectivity and other observables. 
The numerical results for the reflectivity |rp (cf. Eq. (73)) are presented in Fig.js] 





X [|im] 



A. [|im] 



FIG. 5: Reflectivity \r{io,kx)\'^, cf. Eq. (73), as a function of the wavelength. Left: angle of 



incidence 5°. The other parameters are as in Eqs. (57)-(58) apart of the grating period which 
is 0.87 /xm. Solid lines: exact values, dashed lines: values obtained in the resonance dominance 
approximation. Dotted line (in the left figure and coinciding with the dashed line): values obtained 



by truncation of Eq. (31) to n,m = 0,-1, see the main text. Right: angle of incidence 86° with 
changes in the parameters to e/y = 1 and A = 0.8 fim. Solid line: exact values, dashed lines: values 
obtained in the resonance dominance approximation and shifted by 6.4 nm respectively towards 
smaller wavelengths. 

With the indicated parameters of the grating the resonating guided mode Hk^K has k = —1 
and K = 0. Our results are compared with numerical results obtained by applying transfer 



matrix techniques after truncating the system of equations (31). This is similar to what 
is used in the rigorous coupled wave analysis (RCWA), cf. Ref. [19j. Convergence was 
typically achieved with 20-30 channels - slower than in the TE case. In the left part of 
Fig. |5] we also compare with the results of truncation of the system of Eqs. (31 ) to just two 
channels, n,m = 0, —1. 

As in the TE case we have also examined the structure parameters for which the Fano 
interference between the resonance and the background is particularly pronounced. This is 
shown in the right part of Fig.[5j We observe a good qualitative agreement in both graphs. 
The same quality of agreement has been found for other values of the angle of the incident 
hght. 

At the same time it is seen that on the quantitative level the comparison with the exact 
results for TM modes is less satisfactory than for the TE modes. For example in the 
particular resonances of Fig. [s] the discrepancy in the resonance position is —3.9 nm (left) 
and —6.4 nm (right) as compared to 0.6 nm and 0.7 nm in the TE case. The discrepancy in 
the width in the left figure is 0.639 nm between the exact results and resonance dominance as 
compared to just —0.013 nm in the TE case. What is particularly notable is the discrepancy 
with the two channel truncated results where an almost perfect agreement for the TE case 
has been obtained. 
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We tend to attribute these discrepancies to the pecuhar TM boundary condition at the 
borders of the grating layer. These conditions for the guided mode channel are violated in the 
resonance dominance approximation which we use. We have verified that the discrepancies 
decrease with decreasing contrast of the grating or decreasing size of the grating interval. 

We have also noticed that resonance dominance begins to fail in the TM case for regular 
binary gratings with very small or very large duty cycle for which I73I ~ I72I ~ |7i|. 

A more systematic analysis is required to clarify the above issues. In this respect let 
us note the following. The general framework for the resonance dominance approximation 
is provided by the Feshbach projection operators formalism as it is outlined in Appendix 



VIA There the exact resonating state is one of the eigenstates in the Q subspace. In our 



approach we have approximated this state by neglecting the coupling between the Fourier 



components with n k, c.f.,Eq. (64). Using an expansion of the modes in terms of 
Bloch waves rather than plane waves appears as a promising tool to improve the treatment 
of the TM modes. 



IV. INTERACTION OF FESHBACH RESONANCES 
A. Overlapping TE resonances 

1. Formal development 

In this section we will extend our discussion and address the issue of overlapping res- 
onances and their interaction. For photonic crystal slabs of the structure shown in Fig.[T] 
overlapping resonances always exist at illumination close to normal incidence. This has its 
origin in the time reversal symmetry. It is straightforward to verify that with En{kx, kf,z) 



also E*_^{—kx, —kf,z) solves the system of equations (10) and satisfies the boundary con 



dition (11). This implies that resonances occurring for incident light with k^ always 
involve two guided modes, i.e. a resonance with n = u is always accompanied by the reso- 
nance for n = —V. In the resonance dominance approximation this conclusion follows readily 



from the resonance condition (22). If it is satisfied for a guided mode Er^^iz) and n = v sX 
fca; — )■ it will also be satisfied for the same guided mode and n = —v. The physics of this 
condition is that the incoming light with /cj, — )■ is coupled by the grating to both right 
= kx + i^Kg) and left {(3 = k^ — uKg) propagating guided modes having the same 2;-profile 
{z). 

We consider in detail this class of overlapping resonances, i. e. we include together with 



^■qo (z 



Ey also E^y resulting in a system of equations with 3 components (cf. Eqs. (46), (47)) 



Eo{z) = e'^^\z)-uj^ [ dz'go{z,z'){e^,E,{z') + e,E.,{z')) , (76) 
^(i. ^ rFV2^m{^) I dz'S,,{z'){ei,E,{z') + e±2uE^u{z')). (77) 



E^Jz) 



Here we have made use of the peculiar property of the TE modes that the value of u gives 
rise to shifts of the the eigenvalue rjo) but does not affect the eigenmode Sj^^ (cf. Eq. (16)). 



According to Eq. (77) the resonance waves E±,y{z) differ only in their (resonance enhanced) 



strengths a±,^ (cf. Eq. (48)). Up to these unknown normalizations they are given by the 
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guided mode S, 



vo 



(78) 



As in the case of an isolated resonance (cf. Eq. (79)), the reflection amplitude is determined 



by the asymptotics of Eq{z) in (76). Using Eqs. (13) and (78) it is easily seen that Eq. (76) 



is obtained from Eq. (46) by replacing e_y by e^y + a^y e^. Accordingly the expression 



(54) for the reflection amplitude of an isolated resonance is replaced by 



2k: 



(79) 



It remains to determine the unknown strength parameters a±y. To this end one inserts 



Eq. (76) into the two equations (77) and obtains after using (78) a 2 x 2 linear system system 
of equations for the unknown variables (J±y, 



( (^+u \ 













with the coupling matrix W 
W = 



W 



(80) 



\2 , \^ 12, ,4 ■ 



(81) 



The diagonal elements of the matrix W contain the complex valued self interaction E defined 



in Eq. (|52|). The off-diagonal elements contain a direct coupling via e±2y between the guided 



modes with 



V 



(82) 



as well as the indirect coupling terms e^S and e^^^S via the extended mode. In terms of the 
matrix elements Wij, the inverse of W is given by 



1 



W+W. 

where the eigenvalues of W are given by 



W22 -W12 
-W21 Wii 



(83) 



± J(2A;,z/irg)2 + a;4(^(|,^|2^2S-|e2^|cos(^,y)' + |e2.|2sinV.^')- (84) 



We have introduced the relative phase between e2y and e 

r*2 



^2yel 
I 9 1' 



(85) 



which will be seen to distinguish different types of interacting resonances. The importance of 
the relative phase between the different Fourier coefficients has been previously emphasized 
by Barnes et al. [38j in the similar phenomenon of surface plasmon resonance off metallic 
gratings. Here it is not the relative phase (~ £2/61) which matters. Rather the phase tp^ 
results from the interference of the 2-step process via the extended mode and the one step 
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process connecting directly the 2 guided modes. The combination of the strengths which 
determines the reflection amplitude is easily calculated 

(T,e_, + (T-^e, = 2u;^|e,|'C(+) ^^ , Wo = + kl + v^KI + \e2u\ cos ^,u^V . (86) 

This is the central result of this section. We will use it in the next section to analyze possible 
patterns of overlapping resonances. 



2. Reflectivity of overlapping resonances 



The equations (79) and (86) show that the zeros of the real parts of Wj^ and W_ define 
the positions of the two resonances while their widths are given respectively by Im Wj^ and 
Imiy_. According to Eq. (84) the sum of the resonance widths and the midpoint of their 



positions are not affected by the interaction of the resonances. It follows that broadening of 
one of the resonances is accompanied by narrowing of the other. This is reminiscent of the 
phenomenon of motional narrowing [39] as well as the phenomenon of superradiance [IQ] 
and subradiance [H] for coupled emitters. The distance |Re(VF+ — WJ)\ between the two 
resonances can either shrink or expand depending on the structure of the GWS. If 



Re 



\(i2u\ COS 



t2y\ sm V9, 



V'\ <0, 



(87) 



and if the distance \4:kr^vKg\ between the non-interacting resonances is sufficiently large 
this distance decreases due to the interaction. The condition (87) is actually satisfied for 
the choice of the parameters (58) where |ImS| ^ 2.5|ReS| and = 0. Obviously, for 



sufficiently large values of the term proportional to V repulsion of the resonances results. 
Below we will discuss such a case. For kx = the distance trivially has to increase or to 
remain as a consequence of the interaction. In the regime in between these limits the 



change in distance depends on the details of the various quantities in (87). It is remarkable 
that the system of overlapping resonances can be tuned to exhibit either "level" repulsion 
or attraction by variation of an external parameter (kx)- A peculiarity of the interacting 
resonances is the appearance of a zero in the resonance amplitude, i.e. Wo{u,9) = in 



Eq. (86). The presence of this zero may distort significantly the shape of the refiectivity 



or obscure the presence of two resonances if the distance between the zero of Wq and the 
resonance position of W- is of the same size as or smaller than the corresponding width. 

Finally, the ratio of the kinematical term {2kxi'KgY and the interaction induced 2nd term 
in the square root of Eq. (84) controls the transition from overlapping to isolated resonances. 



For sufficiently large k^ the interaction induced term, i.e. the off-diagonal matrix elements 



of in (81 ) can be neglected and two isolated resonances are described by W . 



We now will give a brief overview of the possible structures of the resonances generated 



by Eq. (86) and we will discuss their dependence on the properties of the GWS. In these 



qualitative studies, but not in the numerical results, we disregard the contribution ro(a;, k^) 
to the reflection amplitude r{u,kx) (cf. Eq. (79)). We first consider the case when \e2u\ = 
for which ip^^ is ill defined. This is realized for a grating with 50% duty cycle. According to 



Eq. (84), for kx ^ Q the transition from two separated resonances to one resonance takes 
place. Ai kx = ^ only one resonance is present with a width twice as large as that of an 
isolated resonance (W- is exactly canceled by Wq), while at arbitrarily small non zero angles 
two peaks will appear due to the small imaginary part of W^ . 



20 



This is illustrated in the left part of Fig. |6] for the case where the resonance conditions are 
satisfied for modes with v = ±1. In the resonance region leipcu^S ^ (0.07 + O.lGi) /xm"^. 
The distance between the resonances is thus smaller than the wide resonance width and 
indeed no clear separation of the resonances is observed until the incidence angle \d\ is 
increased to values of the order of or larger than 

1^1 > -leilVlmS/iTg ^ 0.3°. 



One can also see the effect caused by the zero of Wq in Eq. (86 ). It produces the dip in the 



combined peak in the left part of Fig.|6| The appearance of the zero can also be interpreted 
as arising from a cancelation between the contribution from the two eigenstates of W (81) 
associated with the eigenvalues W±. This destructive interference of the contributions from 
the two eigenstates causes a transparency window within the resonance and is closely related 
to the "EIT" phenomenon (electromagnetically induced transparency), cf., Ref. |12]. Thus 
for an appropriate choice of parameters guided mode resonances provide yet another classical 
analog of EIT, cf., |l3l E]. 




1.43 



1.45 




FIG. 6: Reflectivity as a function of the wavelength for the angle of incidence 9 = 0.05°, structure 
and parameters as specified in Eqs. ( 57p8 ). Solid lines: exact values, dashed hnes: values obtained 
in the resonance dominance approximation. Left: duty cycle d/A = 0.5. Right: duty cycle 
d/A = 0.75. 



We now consider the case where and el are, up to a sign, in phase, i. e., for 



0, TT 



(88) 



Then Eq. (84) simplifies for small k^, 



W. 



± 



t2v\ COS 



^9) 

When the direct interaction V is strong in comparison to the self coupling S one will have 
a clear separation of the two resonances. Since V is real one of the resonances will still have 
a vanishing width at fc^^ = while another will have the full width. This is illustrated in 
the right part of Fig. |6] for the duty cycle of 75%. The shape of the reflectivity is drastically 
changed as compared to the left part of this figure showing two well separated resonances 
with widths which differ by a factor of 350. 

It should be noted that any simple binary grating of arbitrary duty cycle (cf. left part 
of Fig (|T|) are symmetric under x-mirror reflection over a plane defined by x = Xq. Any 
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grating with such symmetry will either have = or will fulfill the condition (88). Since 
the incident field is symmetric under x-mirror refiection symmetry only for fej. = we may 
interpret the vanishing width at normal incidence for such symmetric gratings as a symmetry 
selection rule. Thereby it is only possible to excite the combination of guided modes which 
is symmetric under the x-mirror refiection. Interacting resonances with similar widths at 



small or vanishing angle of incidence can be generated only if the relative phase ip^ (85) 
deviates significantly from or vr. 
We consider the case 



±- 



TT 



(90) 



and realize this value by choosing the profile of the grating shown in Fig. [7] with parameters 
A = 0.81 /im, e^in = l, e„,ax = 4, rf/A = 0.2 . (91) 



FIG. 7: Profile of the dielectric constant in an elementary interval of the grating layer with grating 
period A. 

The structure of the resonance curves depend on the strength a of the direct interaction 
of the 2 resonances in comparison to the indirect interaction via the extended mode 



a 







V 


e 


2 

u 




|E| 



(92) 



With the above choice of the parameters (91 ) the direct couphng dominates, i.e. a 3> 1 , and 
we find for small 



W± = r]^ + kl + u^K^ + eJ V S ± 



^ ^2. + . 21 IT/ 

2uj^\e2u\V 



(93) 



The direct interaction V generates a repulsion of the two resonances; to leading order in a, 
their widths are equal that of a non-interacting resonances. Also in this limit the numerical 
results confirm our analytical analysis as is seen in Fig. [81 



3. Contour plots, hand gaps and curvatures 

Contour plots of the refiectivity in the ^ - A plane offer a convenient way to survey the 
scattering in different kinematical regimes including the separated and strongly interacting 
resonances. In Figs. ^ and ( 10 ) such contour plots are shown which correspond to structures 
and parameters of Figs. ^ and (|8|. Quantities which are important for applications and 
which can be extracted from the contour plots are the band gap, i.e. the difference in 
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FIG. 8: Reflectivity as a function of the wavelength for the angle of incidence 9 = 0.005°, grating 
structure shown in Fig.[7] and parameters (57), (58) and (91). Solid lines: exact values, dashed 
lines: values obtained in the resonance dominance approximation. 



wavelength between the peaks |rp = 1 in the limit — and the curvature of the |rp = 1 
curves at these points. We will calculate these quantities and establish their dependence on 
the parameters of the GWS. 

An analytic understanding is possible only under the assumption that the coupling of 
extended and guided mode is weak, i.e., [eyC*^"*")] a; <^ 1. This condition is satisfied in all of 
the examples to be discussed. It implies that values of |r(ci;, kx)\'^ of the order of 1 can be 
realized only if uj'^\Wo/W+W_\ 1, i.e. if Reiy+ or ReW_ ^ 0. In the following analytical 
studies we therefore replace the contour lines of the reflectivity by the contour lines of the 
equations 

ReW±{uj,e) ^0. (94) 



Fig. 10 displays the level repulsion arising from the direct coupling term in Eq. (93). We 



know from Fig. [6] that the strength of the level repulsion is correctly reproduced in resonance 
dominance. To calculate band gap and curvature we drop the sub leading terms proportional 
to S in Eq. (93 ), expand around uq = 27r/Ao, the zero of ri{ijJo)+h'^Kg and obtain the following 
expressions for the band gap 6X and curvature k 



5X 



-V'o 



dX 



An 



1 ± 



;/2 \2 



(95) 



22 nm, and for the curvature of 



These estimates yield for the structure of Fig. [10] SX 
the 2 contour lines k ~ 0.3 ± 6.8 nm , and agree with the numerical results of Fig. 10 



Analytical results for the contour plot on the left hand side of Fig. |9] are harder to obtain 
due to the presence of the zero (cf. Fig. [6]) of the resonance contribution to the reflection 
amplitude (cf. Eq. (86)) for the parameters (58) of the GWS. We note, that €2 = for the 
choice of the duty cycle d/A = 1/2. As discussed above, the transition from isolated to 
interacting resonances takes place if 



2u9K„ 



1. 



(96) 



With leipw^S ^ (0.06+i0.16)/im ^ this condition is satisfied for 6 = 0.15° . For larger values 
of 6 we expect a linear dependence of the wavelength on the angle of the incident light in 
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FIG. 9: Contour lines of the reflectivity in the 9-\ plane. Structure and parameters are specified 



in Eqs. (57 58). In the right part of the figure the duty cycle has been changed from 50 % to 75% 
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FIG. 10: Contour lines of the reflectivity in the 9- \ plane. Structure and parameters are specified 



in Eqs. (57) ,(58) and (91). 



agreement with the numerical results. For significantly smaller values of 9 we consider the 
resonance associated with W+ which is not directly affected by the presence of the zero of 
Wq. We identify the lower branch of the contour lines of the reflectivity with the contour 
lines Re Wj^ = 0. Expanding as above W+ around uq we obtain 



2Kl 



110 nm . 



(97) 



i.e., these estimates account for the order of magnitude difference of the curvatures in Fig. 10 
and the left part of Fig.|9| 



4-. Electric fields of overlapping resonances 



So far, we have concentrated our discussion on one observable, the reflectivity. An inde- 
pendent observable is the resonating electric fleld with components E±u (cf. Eq.(78)). With 
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only minor modifications the above analysis of the reflectivity can be applied. In terms of 
a±, the guided mode part of the electric field is given by (cf. Eq. ) 



EJx, z) + E^Jx, z) 



-i ^e_|_ cos vKgX + e_ sin uKgX^ e^^^^^Er^Q {z) 



with 



(98) 



(99) 



The calculation of the two components via Eqs. (80) and (83) is simplified if we choose to 



be real (which always can be achieved with an appropriate choice of the origin of x in the 
integral ([6])). It follows immediately from Eqs. (86) and (79) that e+ is proportional to the 
resonating part of the reflection amplitude 



2k^ {r{uj,k^) 



'■CM 



(100) 



This expression is valid for any value of 6'. It applies also to the the case of isolated reso- 
nances. Up to the normalization, the coefficient e+ is given by the (resonating part of the) 



reflection amplitude. As seen in Fig.[TTJ e+ vanishes at uq and exhibits the same asymmetry 

The absolute value of e_i_ 



around the corresponding wavelength as the reflectivity in Fig.[6 



is related to the half width of the isolated resonances (cf . Eq. 53 ) by 

|e+| = 2a/ k-/t |r(w, k^) - ro{uj, k^)] 



(101) 



i^e. the strength of this component of the electric field increases with decreasing half width 
f . In region III of the GWS (Fig.[l|, the strength of the electric field proportional to e+ can 
be larger than the incident field by by up to a factor of 10 (cf. Fig.|3]). The e_-component 



(102) 



e_ = ^ ^ — ( - ^k^uKg + z|e2j.|a;Vsinv9e) 

of the resonating electric field differs from e+ by the factor W^^ . For 

V9 = 0, vr 



the appearance of Wq in Eq. (86) is necessary to keep the reflection amplitude finite by 



canceling the zero of W- (cf. Eq. (89)). In turn this implies that e_ is singular for ^ — )■ 
and uj = ujq where ujq denotes zero of Wq 



W^M = rioiul) + v'^Kl + ujl\e2u\V{ujl) = 0. 



(103) 



For values of {9, u^) sufficiently close to the singular point (0, u^) in the 6',a;^ plane, e_ is 
given by 



(104) 



Considered as a function of u, for fixed 9, the component e_ has a Lorenzian (Breit Wigner) 
shape. It reaches its maximal value a.t u = uq — 6'^ReS_ 



4:x Tk- 



IS-I 



MuoKq 9 \ ImE_ 



(105) 
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FIG. 11: Electric fields e+ (left) and e_ (right) (cf. Eqs. (100) and ( 102| )) in units of the amplitude 
of the incident field as a function of the wavelength for the interacting resonances of the left hand 
side of Fig.joj Solid lines: real part, dashed lines: imaginary part, red (color online) 9 = 0.05°, blue 
(color online) 9 = 0.01°. The values for e- at 9 = 0.01° have been divided by a factor of 5. 



As Fig. 11 demonstrates the increase in the electric field strength with decreasing angle 9 
of the incident light is accompanied by a decrease of the half width 6'^ImS_ of the Breit- 



Wigner function (104). In the homogeneous layer III (cf., Figs.Hl^ the resonating part of 



the electric field reaches values which are up to a factor 3/6* [°] larger than the field of the 
incident light. 



B. Overlapping TM resonances 

1. Formal development 

We now consider interacting resonances in the TM case. We follow the developments in 



Section IV and consider the case when there are two eigenvalues r]i = rjn^ y^ and rj2 = f]n2,u2 
with values close to u"^. As the comparison of Eqs. (16) and (38) shows, unlike in the TE 



case, the guided modes 'Hn,^, depend on the index n which makes the following analysis 
more involved. The same complications will occur in the TE case if overlapping resonances 
associated with different eigenvalues rjn^ ^ considered. 



We retain both resonating terms in the system of equations (40) and write coupled equa- 



tions for the resonating and the zero components. We will then use the resonance dominance 
approximation retaining only a single term with the guided mode l-ini,Vi{z) in the sums over 
iinEq.(|42l) 



Hn2{z) 



HI-'\z)+ [ dz'g{z,z'){eo,n,HnAz') + eO,n2Hn2{z')) , (106) 

[ rf^'H„,,,,(z')(0n„O^o(^') + 0ni,n.^n.(^')) , (107) 

w - Vl Jig 

[ dz'nn2,U2{z){Qn2flHo{z) + Qn2,rnHnAz')) ■ (108) 

Jia 



'^n2,U2 \Z) 



As in the case of the single resonance the components if^j and iJ„2 are proportional to the 
guided modes 



Hnii^Z^ 0'niT~^ni,ui(^Z^ , ^Tyjg (^) '^112'^^ 



(109) 
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but the equations determining o"„.'s show mixing of the two resonating modes. Inserting 



(109) into (107) and (108) together with (106) we obtain closed equations for the field 



enhancement coefficients 



W ( ^"^ 



-7 C^^^ 



(110) 



where, in analogy with the TE case (Eqs. (80), (81), we have introduced the coupling matrix 



W 



'1n2-ni^ 7-rii7n.2^«2,ni 



"2 



and 



111) 



(112) 



We have used the replacement explained in Appendix VI C — ^ In-mQnm with the 
notation _ ^ 

Oram (^z 

d, + {K + nKg){K + mKg). (113) 

The quantities 



(114) 



dz' I dz"'Hn,,uXz')<dn„o9o{z',z")Qo^n,T-Ln,,u,{z") , Z,J = 1,2. (115) 



are generalizations of the definitions ( 71 ) and ( 72 ) for isolated resonances. Here a new 
quantity appears 

V = dz''Hni,uiiz')Qni,n2'Hn2,U2iz') (116) 
J I, 

which controls the direct coupling between the two guided modes. An indirect coupling 
between the guided modes via the interaction with the extended mode is generated by 
Sni,n2 = Sn2,ni- The quantities E^.^^ii are the self-interactions of each guided mode via the 
extended mode. Note that for the dielectric structure in which = = but with 
7ni-n2 7^ we obtain the standard avoided level crossing problem, i.e. the non radiating 
photonic band gap case. 

The eigenvalues of W and the field enhancement coefficients are given by 

± = uj^-(^,+^^)/2 (117) 



W. 



±A/(/il-/"2)V4+(7 ni-n2 

v + )(7n2-niV' + 7— rti7n2'^n2,"iy ' 



7ni(w^ - /i2)Ci|^ + 7n2(7ni-n2'^^ + 7ni7-n2 ^ni ,n2 )Cnt^ 



' 711 



w.w. 



a. 



7ni(7 n2-ni v + 

7-ni7n2^n2,ni )d|^+7n2(w'-/il)C 



"2 



"■2 



(118) 
(119) 



Using Eq. ( 106 ) and the asymptotics of the Green's function (cf. Appendix VI B ) the 

^ {^n.l-uAV + «^n27-n2Ci:)) (120) 



reflection amplitude is given by 
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2. Overlapping resonances close to normal incidence. 



To proceed with the analysis we now specify to the special case of the vicinity of = 0. 
We note that a solution of Eq. (39) with a given fc^, and n is also a solution with k^ — )■ —k^; 
and n — )■ —n. This means that ai k^ = the resonance condition (43) will always be satisfied 
by a pair of guided modes. Let us consider such a resonating pair with the eigenvalues 
Vn,u = V-n,u- For such a pair ni = n = — 77-2 and at kx = all the components of C^~^^ and of 
S are identical and given by 

c(+) = cw= / dz'[d^mnAmd.'Ht\z')] , 
Jig 

S = Sn„n, = / dz'dz''[d,mnA^%d,'d,,,go{z,z'')][d,.HnA^'')], (121) 

Jig 

V = [ dz' {[d^'-HU^')]' - {nK^n-HnA^')]'} . 

Jig 

For small deviations from fc^; = we find 

/oo 
dz^o{z)nlAz) , (122) 
-oo 

where we have used Eq. (38) and 

dz 'Hn,v~^l T^n^u ■ (123) 



dkx J—r^ ' dk 



Identifying for simplicity S and V with their values at k^ = 0, we approximate Eq. (117) 
by 

W± = [lo^ - r]n,u - bn\^^ (124) 



±y {2CkxnKgy + {{-fnl'^i:'^ + \72n\cosipyV) +|72nPsin ip^V^ 

with ^ 

e'^'< = . (125) 

l72n7nl 

The reflection amplitude is given by 

r(cu, fc^) = ro(a;, fcx) ^ , VV o = ^ |72n I cos v^-yV/ . (126) 

The above expressions for W± and r(a;, k^) have the same form as in the TE case (cf. Eqs. 
(84), (86)). We again find that the interference of the contributions of the two eigenvalues 
W± of the matrix W gives rise to a zero in the resonance part of the reflection amplitude. 
Also the role played by the phase y?^ is identical to that of the TE phase In particular 
for values 

= 0, TT 

the width of one of the resonances (corresponding to WA) goes to zero at kx = 0. This is 



illustrated in Fig. 12 In Figs. [13] we show an example of resonance dominance compared 
to the exact results. 
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FIG. 12: Contour lines (exact results) of the reflectivity for TM polarization in the 6-X plane. 
Structure and parameters as given in Eqs. ( 57p8 ) apart of the grating period O.SJ/jm and 50% 
duty cycle (left figure), 75% duty cycle (right figure). 




FIG. 13: Reflectivity for TM polarization as a function of the wavelength at 9 = 0.05° for the 
same structure and parameters as in Fig. 12 with 50% duty cycle (left figure), 75% duty cycle 



(right figure). Solid lines: exact values (vertical cut in Fig. 12), dashed lines: resonance dominance 
approximation (cf., Eq.(125) with ip^ = 0). 



Note that the dependence on the duty cycle of the TM results is qualitatively very 
different from the TE case, Figs. [6]and[9j The apparent band gap at = closes at a very 
different value of the duty cycle. The source of this difference is clearly seen in our theory. 
Let us define the band gap as the difference between the two resonant frequencies at normal 
incidence. This corresponds to the difference between the real parts of and W^. When 
the band gap size is significantly larger than the FWHM of the wider of the two resonances, 
the two resonances are well separated and thus the band gap becomes clearly visibly as in 
the left part of Fig. [T2j For the band gap of the order or smaller than the wider resonance 
width the resonances overlap and appear merged. 

According to Eq. (84) and (125) the difference between W+ and W- for = 0, vr is given 

by 



|VS + |e2.|^) and 2 {\^^\^^ 



1 72, 



(127) 
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for TE and TM respectively. When the duty cycle is close to 50% we have vanishing \e2u\ 
so that the direct coupling term V does not contribute in the TE case. At the same time 
72„ does not vanish in the inverse Fourier transform matrix approach which we adopted 
as explained in Section |III B[ This difference between the TE and TM is however not 



as 



important as the difference in the expressions for E. They are respectively given by Eqs. 



(52) and (121) with the reference to expressions (149) and (169). Comparing we notice a 



peculiar term present only in the TM case. This is the last term in (169). It is proportional 
to the 5-function and produces a "contact" term 



- 1/70 / dz'[d,,nnAz'f 



(128) 



in the expression for S in the TM case. We call this term "contact" since in it the mode 
functions dzT-Lnviz) "interact" at one point. This is in contrast to expression (84) for the 



TE or (125) with only the first term of (169) for TM. In those cases the mode functions 



interact via a finite range interaction caused by the propagation in the extended mode. 

The contact term is real and therefore contributes to the band gap. Numerical evidence 
shows that this is the dominant contribution for the 50% duty cycle and is significantly 
larger than the wider resonance width given by the imaginary part of S. This is the reason 
the apparent band gap is present at this duty cycle in the TM case. In order to complete 
the argument and to understand why the real part of E does not give rise to the same effect, 
we also compare the corresponding quantities for TE and TM case. To estimate the above 
contribution to the TM ReS we write in the thin grating case 



"9 ' '-n,u 



We compare this estimate with the corresponding TE result 



leiyPw^ReS 



6' 



recalling also that we have to divide the TE result by drjQ/duj'^ ~ 3.5 before we can interpret 
is as a shift. 

As one changes to higher or lower values of the duty cycle the second term, the direct 
interaction V begins to play a role in Eq. (127). Here again the explicit expressions of V, 
cf., (82) and (121) show the clear difference between the TE and the TM waves. In the 
former V is real, positive and therefore contributes to the increase of the band gap. For the 
TM case is a difference of two positive terms and can be either positive, i.e. increasing 
the gap or negative, i.e. leading to its decrease. In fact we find that the apparent TM band 
gap closes at 75% duty cycle, cf. the right part of Fig. 12 as well as at 17 %, unlike in the 



TE case where the band gap closes only for a single value of the duty cycle. 

The overall agreement between the resonance dominance and the exact results in the 
TM case is not as satisfactory as for the TE case. As an example we show in Fi g. [Til the 
TM results for the same complex grating as in the TE case, cf.. Fig. g and Eq. (|9l) with 



\g = 0.87/im. It is seen that the band gap is significantly smaller and the relative widths 
of the two resonances are different in the resonance dominance approximation as compared 
to the exact results. In our understanding there is a number of effects in the TM case 
contributing to such a discrepancy vis-a-vis an excellent agreement which is obtained in the 
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TE case. Firstly, it is known, cf., Ref. |33] that numerical convergence as a function of 
number of truncation orders is much slower in the TM case than the TE case. We note that 



the grating parameters as in Eq. (91) lead to I72I ^ |7i| raising the relative importance of 
the higher order components. 

Secondly, we believe that the issue of the special boundary conditions (cf., the end of 
Section IIIB) in the TM case is not accounted properly in our way of implementing the 



resonance dominance approximation. We plan to address this and related problems of the 
TM polarized scattering in the future work. 
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FIG. 14: Left: contour lines (exact results) of the reflectivity for TM polarization in the 6-X. 
Right: reflectivity for TM polarization as a function of the wavelength at 6 = 0.005°. Solid line: 
exact values (vertical cut in the left figure), dashed line: resonance dominance approximation (cf., 



Eq.(125). Same structure and parameters as in Fig. 12 apart the grating parameters which are 
given by Eq. (91) with = 0.87/xm. 



C. Comparison with coupled resonances in quantum mechanics 

The physics of coupled guided mode resonances of light scattering off photonic crystal 
slabs discussed above is closely related to the physics of interacting (overlapping) resonances 
or bound states in quantum mechanics. This connection is obvious in the limit of vanishing 
coupling to the extended mode. Then the matrix W (81) is hermitian and the two (real) 
eigenvalues W± of W can be written as 

W± = ^ [W, + W2± ^/{W,-W2y + Au^\e2u\'V^] , 1^1,2 = Vo + {K ± JyKg^ . (129) 

It's counterpart is the quantum mechanical two level system defined by the Hamiltonian 
matrix 

with the energy eigenvalues E± 

E± = \ (^1 + ^2 ± ^{Ei - E2Y + 4t;2 j . (131) 
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The level repulsion, i.e., increase of the gap of the two eigenvalues with increasing or 
is the characteristic property of both systems. To establish the correspondence in the 



general case is more involved since the matrix W (81) is not hermitian. The non-hermitian 
contribution arises due to the non- vanishing imaginary part of S. The (identical) contribu- 
tions of Im(E) to the diagonal elements of W is due to the decay into the extended mode. 
Similarly the imaginary part of the off diagonal elements, i.e., in the coupling between the 
m = ±z/ modes, S appears since this coupling is generated by intermediate excitation of the 
extended mode. Nevertheless W is not a general complex matrix. The conservation of flux 
restricts the non-hermitian part. Separating hermitian and non-hermitian contributions 

W = Wh- iWnh, 

one easily verifies that the non-hermitian part possesses the following structure 

Wnh = w^'w, 

with 

w = w^VlmS I 

In quantum mechanical scattering on many body systems such as nuclei or atoms, the 
effective Hamiltonian whose eigenvalues are given by the resonance positions and widths 
exhibits exactly this structure of its non-hermitian part [3], [IS], [IS]- 



V. CONCLUSION 

The focus of our studies has been on the nature and the properties of resonances which 
are generated in scattering of light off photonic crystal slabs. We have established that 
these resonances are Feshbach resonances, i.e., they are the optical analog of a particular 
type of resonances formed in quantum mechanical scattering of particles off atoms or atomic 
nuclei. In quantum mechanics, bound states of many body systems in the absence of the 
coupling to the projectile are the progenitors of the resonances. In classical optics the role 
of the progenitors is played by the guided modes of the dielectric slabs in the absence of 
the coupling to the extended modes where the grating region is replaced by a homogeneous 
layer with an effective dielectric constant. Turning on the respective couphngs, bound states 
and guided modes are converted into resonances. Close to the resonance, the interaction 
between projectile and target is dominated by resonance formation. We have shown in 
our investigations that resonance dominance is not only an important concept but also a 
quantitative tool in the analysis of scattering of light off photonic crystal slabs. 

In the formal part of our studies we have employed techniques developed in the context of 
nuclear reactions to establish the resonances observed in scattering of light off photonic crys- 
tal slabs as Feshbach resonances. Despite of the common framework of resonance dominance 
we have identified significant differences in the implementation for TE and TM polarizations. 
We have derived expressions for the observables, the reflection and transmission amplitudes 
and the strengths of the magnetic and electric fields. The various shapes of the reflectiv- 
ity as function of the wave length are properly reproduced. Also strong distortions of the 
Breit-Wigner form are correctly described in resonance dominance and shown to result from 
the interference between background scattering and resonance formation. In comparison 
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with results of exact numerical evaluations the accuracy of resonance dominance has been 
established to range, depending on the particular application, from a few to 15 % . 

In comparison with exact numerical methods, the power of the resonance dominance ap- 
proach rests upon the analytical formulation of the observables. It offers the possibility to 
study in detail the dependence of reflection and transmission amplitudes and the strength of 
the electromagnetic fields in terms of the properties of photonic crystal slabs and the kine- 
matics of the incident light. We have demonstrated this potential of resonance dominance 
in the analysis of the rather involved phenomena appearing in the excitation of interact- 
ing resonances. In particular, we have been able to identify the source for the significant 
differences in the interaction of TE and TM resonances and the relevant parameters of the 
grating layer which control the shape of the reflectivity of overlapping resonances. We have 
derived analytical expressions for the band gap and the related curvature as a function of 
the angle of incidence and have study analytically the weird behavior of the electromagnetic 
field enhanced by the overlapping resonances. Needless to say that for design issues this 
analytical procedure can be reversed, i.e. the properties of the photonic crystal slabs can be 
identified which optimize certain requirements on observables. 



VI. APPENDICES 

A. Feshbach projection operators for light scattering 

Since Feshbach's seminal work |lH3] on the theory of nuclear reactions, the main tool for 
such investigations is the formulation in terms of projection operators acting in an appro- 
priately defined Hilbert space. We now will sketch a reformulation of this approach in the 
context of differential equations describing the classical scattering of light, [H Uni IHj. 

We will first consider the TE case. It is convenient to begin with the coupled equations 



(10). We will rewrite them in the form 

oo 

nm ^m~n 

{z))]Em{z) = u En{z) (132) 

m=— oo 

or in the matrix-operator notation 

hip = u'^ip (133) 

where 

hnm = [-dl + {K + nKg){k^ + mKg)]Snm + (^'^{Snm " em-niz))] (134) 

and 

iIj{z) = {..., S_™(z), E_i{z), Eo{z), E,{z), E^{z), ...} (135) 
We define projection operators acting on the array ip{z) 

Pip = ipp = array obtained from ip by setting to zero all components apart of ipo , 
Qip = ipQ = array with the same components as ip but with ipQ = (136) 

Clearly 

P^ = P , Q^ = Q ,P + Q = 1 (137) 



33 



In a way which is standard in the Feshbach formahsm [T],[3] we can act with P and with Q 
on our basic equation hip = u'^ip. Using the above properties of P and Q we obtain coupled 
equations of the Feshbach formahsm 

hppipp + hpQil)Q = u^ipp (138) 

where we defined 

PHP = hpp , QHQ = hQQ , PhQ = hpQ , QhP = hQp. 

Clearly hpp = /iqo; ^qq is the matrix hnm with n = row and m = column excluded. The 
operators hpg and hqp are matrices consisting respectively of just n = row and m = 
column both without the m = n = element. We can formally solve the second equation 

ifjQ = {uj^ - hQQy^hgpifjp (139) 

and insert into the first equation Converting the resulting equation into an integral form we 
obtain 



Mz) = 4''\^) + j dz'go{z,z'){z'\hpQ\iPQ) = (140) 
= i^o'\z)+ / dz'go{z,z'){z'\hpQ{u^ - hQQ)~^hQp\ijp) 



where tpj^^ is the properly chosen solution of {oj'^ — hpP)il)\^'' = and go{z, z') is the (properly 
chosen) Green's function of this equation, i.e. (w^ — hpP)go = 1. The meaning of "properly 



chosen" for the present problems is explained in Appendix VI B We have also assumed that 
there is no unperturbed solution in the Q sector, i.e. ipQ = ^ for hpq = 0. 

So far we made no approximations. Let us consider the eigenfunctions of hgg 

hQQ<pu = riu4>u (141) 

and assume that the geometrical parameters of the optical system and the kinematics of the 
light scattering are such that there exist a range of discrete values of fji, (guiding modes) . We 
will further assume that hpp has a range of continuum eigenvalues (radiating modes) which 
overlap with the discrete modes of hgg and moreover that lies within this combined 
range. In these circumstances the "off diagonal" terms hpq and hgp couple the guiding 
modes to the radiating modes turning the former into Feshbach resonances. Note that we 
use fjiy to distinguish from rji, used in the TE section of the main text. 

A generalized version of the resonance dominance approximation used in the present 
work amounts to using one or few resonant terms in the full spectral decomposition Green's 
function in the Q sector 

^V = E¥^+G. (142) 

where Gc is the contribution due to the continuum eigenfunctions. Using as an example two 
resonating eigenstates we obtain 



1 ^, |0l)(0l| , |02)(02 



(143) 
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which means (cf., (|139|) that 

IpQ = 0-101 + Cr202 



{u?-r]i) ^(0i|/iQp|?/^p) , z = 1,2 



Following the standard route one inserts this into ( 140 ) 



%ljp{z) 



^o^H^) + dz'go{z,z'){z'\hpQ\ 

1=1,2 



(144) 



(145) 



then multiplies and integrates by / dz{(t)j\hQp\z) ... with j = 1 and j = 2 and obtains coupled 
equations for cXi 



- r^i - 

"^21 



Sll — Si2 

U"^ -fj2- S22 



0-2 



Here 



= ((/),|/lQp|4+)) , 2 = 1,2 



(146) 



(147) 



In the TM scattering one can follow exactly the same route by simply noting that the basic 
set of equations (31) is already in the form hip = u'^ip with hnm = ©r 



and 



ij{z) = {..., H^U^), H^,{z), Ho{z), H,{z) 



Hm{z), ...} 



(148) 



It is important to recognize that equations ( 146[ ) are based on assumed exact eigenfunc- 
tions 01 and 02 of Hqq, Eq. (141). In reality hqq is almost as complicated as the full h in 
Eq. (134) and further approximations are required. In our treatment of overlapping reso- 



nances leading to (80) and (110) we have drastically truncated the problem and replaced 
the exact 0i and 02 by the eigenfunctions of hu and /i22- In this way the coupling between 
these components via hu and /121 as well as the coupling to other components are neglected. 



In fact we observe that equation (146) is similar in form to the equations (80) and (110) 



apart from the absence of the term V in the off diagonal elements of the matrix W. This 



term approximately accounts in Eqs. (80) and (110) for the direct coupling between the 
approximate guided modes 0i and 02. 



B. Background fields and Green's function 

In this section we will derive the necessary expressions for the background field and the 
associated Green's function. As in the main section we assume that the dielectric function 
eo{z) approaches 1 for 2; —t- —00. 

1. TE waves 



denote by eI^\z) two independent solutions of the homogeneous equation (14). Using 
these solutions the Green's function can be expressed as 



We start with Eq.(12) for the Green's function g^^z^z') in the TE scattering. Let us 



1 



g,{z, z') = --{e{z' - z)E\^\z') E'-\z) + e{z - z')E\^\z) Et\z')) 



(149) 
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where W is the (2- independent) Wronskian 



W='^Et\.)-Ei-\.)'-^. (150) 



This is easily verified by observing that by construction the Green's function satisfies Eq.(12) 



for z ^ z'. Furthermore the pref actor —1/W is found by integrating (12) over z in an 
infinitesimal interval around z'. 

The choice of the two hnearly independent solutions eI'^\z) is dictated by the type of 
the Green function one needs. It is convenient to require that go{z, z') as a function of z for 
the values of z' inside the grating behaves as a reflected (transmitted) wave in the z —00 
{z — )■ 00) region far away from the grating. This leads to the boundary conditions 

lim E^o^\z) = 4{u)e"'^' , lim E^~\z) = toi^)e~'''^" (151) 

2— >— 00 

with kf denoting the z-component of the incident wave vector in the two asymptotic regions 

= wa/ Civ — sin^ 6 , k~ = cja/ ej — sin^ 6 . (152) 

These equations describe electric fields incident from above (Eq) or below (Eq) (cf. Fig.jl]). 
Their normalization is fixed by the amplitude of the incident wave which we choose as 



lim E'o^'iz) = e'^^ ' + r+(u;)e-*'=^ ' , lim (2) = e"*'^^ ' + r^iuj^^ ' . (153) 
Evaluating W in region IV yields 

W = -2iktt+{io) (154) 
Equally well the Wronskian can be calculated in region I yielding the relation 

K{u^)\" + ^Hiu;)\' = l. (155) 



The two solutions Eq'^\z) are related to each other 

Et\^) = 7Jr-M''^''i^) - rrEi;\z)) . (156) 

tg [UJ) 

The Green's function at large negative z is needed in the computation of the reflection 
amplitude 

lim g,{z, z') = ^e-^'^^ Ei;\z') . (157) 
We also need the spectral representation of the Greens-function 

9o{h,z,z)--J dk^ ^2 _ k,2 + • (158) 

Here we have made explicit the dependence of the fields and the Green's function on the 
incident wave vector. Again it is obvious that for z ^ z' the Green's function satisfies 
Eq.(12). With the help of the completeness relation 

/*;B-'K..)^r'-(*:./)^2,.(.-.'). (159) 
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the correct form of the singularity in Eq.(12) is confirmed. Using the decomposition of the 



denominator into principal value and 6 function contributions 

1 P 



the imaginary part of the Green's function is obtained 



eI;\K, z) E^^^\K, z') + z) E^^\K, z' 



Mi 



(160) 



2. TM waves 



The Green's function go{z,z') for the TM case, Eq. (37), is found in the same way as for 
the TE scattering. We denote two independent solutions of the homogeneous equation (i.e. 
of Eq. (37) with the right hand side set to zero) by h'^\z). Then 



9o[.z,z 



Hi'-\z)Hl\z')9{z - z') + H^-'\z')Hl\z)e{z' - 



r(-)r.' 



(161) 



with a z-independent constant v which is found by integrating (37) over z in an infinitesimal 



interval around z'. The only subtlety relative to the TE case is the presence of the dz^{z)dz 



in the defining equation (37). This leads to two modifications relative to the TE case. On 



one hand the expression for v is now 



i^{z)w{zy 



(162) 



where W{z) = Hq^^OzHq ^ — Hq ^dzH^^^ is the Wronskian of the two solutions. On the 
other hand the Wronskian of any two solutions like -ffg^'* {z) is not a constant but obeys the 
so called Abel formula written in our case as '~fo{z)W{z) = constant. Together we have that 
u is indeed independent of z and can be calculated in a convenient point. 

We again choose the solutions Hq'^^z) and h'^\z) such that go{z,z') behaves 
fiected (transmitted) wave when z' is finite and z — ?■ — oo {z — )■ oo). This leads to the 
conditions 



z — 7- — OO : 



H, 



z ^ oo 



H, 



with 
z — 7- oo : 



z) e 



ikz z _|_ ^-\- ^ikz z 



z — > — oo 



H, 



(+), 



z -> e 



+ rn e 



-ik, z 



At z ^ — oo we have 'y{z — j- — oo) = 1/e/ and W{z — )■ — oo) = —2ik^ t^ so that 

^ ~ _ 

2k ^ 



(163) 



. (164) 



(165) 



Using in Eq. (69) the asymptotic expressions for h'^\z) and the Green's function (161) we 
obtain the result (73). 
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C. Dealing with singular expressions in the description of TM resonances 



The derivative terms in the operators B^o and Bqa: contain combinations dztk{z)dz which 
produce (5-function hke singularities when acting on H^{z), 'Hk,K and g^^z^z'). Here we 
will explain how we have regularized such expressions. Let us consider as an example the 



expression (72) for E. Our strategy is to start with narrow but continuous transition layers 
between the dielectric layers. We then transform this expression in such a way that in the 
limit of sharp boundaries i.e. zero transition layers width it will not contain ambiguous 



singularities. Let us demonstrate this for the term of (72) which contains the differential 
operator parts of both 6^0 and 6^0 



dzdz'Hk,K{z)dz-fk{z)d,go{z, z')dz'-f^k{z')dz'Hk,K{z') = 

rf2;(i2;'[a^Hfc,«,(2;)]7fe(z) [9^(9^/5-0(2;, 2;')]7-fc(2;')[^^''^fcA(^')] ^ (166) 
^ |7fcr / dzdz'[d,nkA^)][dA'9o{z,z')][d,mkA^')] (167) 

where we started with infinite integration range assuming that 'y±k{z) tend to zero fast but 
continuously outside the grating layer. In integrations by parts the surface terms vanish 
while the resulting integral has only finite discontinuities in the limit of sharp boundaries. 
Such discontinuities do not cause ambiguities so that the limit can safely be taken. The other 



two terms containing differential operators in ( 72 ) can be regularized in the same manner. 



The above procedure can be conveniently summarized by stating that expressions like 



(72) should be used with the operators Qnm replaced by 

©nm 7n-m[£ S + (^z + nKg){k^ + mKg)] (168) 

where the arrow above the derivative indicates that it acts on the function to the right or 
to the left of it depending on the direction of the arrow. 
We also note a useful identity 

d,d,>go{z, z') = u[dM'"\z)dz'Hl,-\z')eiz - z') + d,,Hi^\z')d,Hl^\z)9{z' - z)] - 

-lhQ{z)5{z - z') (169) 



which is obtained by differentiating Eq. (161) and using (162) 
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